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Abstract 

Introduction: Early rheumatoid arthritis (RA) patients may show rapid radiographic progression (RRP) despite rapid 
initiation of synthetic disease-modifying anti-rheumatic drugs (DMARDs). The present study aimed to develop a 
matrix to predict risk of RRP despite early DMARD initiation in real life settings. 

Methods: The ESPOIR cohort included 813 patients from the community with early arthritis for < 6 months; 370 
patients had early RA and had received methotrexate or leflunomide during the first year of follow-up. RRP was 
defined as an increase in the van der Heijde-modified Sharp score (vSHS) > 5 points at 1 year. Determinants of RRP 
were examined first by bivariate analysis, then multivariate stepwise logistic regression analysis. A visual matrix 
model was then developed to predict RRP in terms of patient baseline characteristics. 

Results: We analyzed data for 370 patients. The mean Disease Activity Score in 28 joints was 5.4 ± 1.2, 18.1% of 
patients had typical RA erosion on radiographs and 86.4% satisfied the 2010 criteria of the American College of 
Rheumatology/European League Against Rheumatism. During the first year, mean change in vSHS was 1.6 + 5.5, 
and 41 patients (11.1%) showed RRP. A multivariate logistic regression model enabled the development of a matrix 
predicting RRP in terms of baseline swollen joint count, C-reactive protein level, anti-citrullinated peptide 
antibodies status, and erosions seen on radiography for patients with early RA who received DMARDs. 

Conclusions: The ESPOIR matrix may be a useful clinical practice tool to identify patients with early RA at high risk 
of RRP despite early DMARD initiation. 



possible, disease remission [3,6-12]. Implementation of 
these recommendations in practice has led to better clini- 
cal outcomes [6,8,13]. 

The choice of the first DMARD has been the topic of 
many trials and guidelines. Methotrexate (MTX) has 
been recommended as the anchor drug because it allows 
for step-up strategies, that is, the addition of other syn- 
thetic or biological DMARDs if adequate response is not 
achieved with MTX monotherapy [3,9,14]. Leflunomide 
is the alternate choice because both drugs seem to have 
equivalent symptomatic and structural efficacy [15]. A 
number of trials have assessed the potential benefit 
of biologic DMARDs as first-line treatment for early 
RA [16-18]. These intensive options have been consid- 
ered more efficacious than MTX in trials with a static 



Introduction 

The care of rheumatoid arthritis (RA) has profoundly 
evolved during the last decade because of new drug thera- 
pies and the early treatment paradigm. RA requires rapid 
referral to a rheumatologist [1-3] and early initiation of 
disease modifying anti-rheumatic drugs (DMARDs) to 
prevent long-term disease consequences, such as diabetes 
or hypertension [4,5]. At the same time, regular adaptation 
of DMARDs depending on disease activity - that is, RA 
tight control - has become an additional standard in RA 
management to achieve at least low disease activity and, if 
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therapeutic strategy during the first year [18]. However, 
in trials adopting dynamic step-up strategies, the overall 
benefit of biologies as a first-line agent remained ques- 
tionable [6,11,12,19-21]. In addition, several economic 
evaluations reported that incremental cost-effectiveness 
ratios of biologies as first-line treatment for early RA 
are usually high and largely overtake the generally 
accepted thresholds [22,23]. These conclusions rein- 
forced the position of MTX (or leflunomide) as the opti- 
mal first-line agent for early RA. 

Despite current recommendations, for some patients, 
MTX may be suboptimal as first-line therapy. Several trials 
have shown substantial structural progression even with 
MTX started rapidly after disease onset. This situation has 
led to the development of the concept of rapid radio- 
graphic progression (RRP), defined as structural damage 
progression of at least five points of the van der Heijde- 
modified Sharp score (vSHS); the cut-off of five points cor- 
responds to the destruction of one small joint and to the 
usually reported smallest detectable difference (SDD) 
[24-26]. The rationale for this threshold has been con- 
firmed in two different studies. In the BeSt trial, patients 
with RRP during the first year of follow-up showed 
increased functional limitations and structural damage 
progression over eight years of follow-up, despite a tight 
control-based therapeutic strategy [27]. These results are 
consistent with another study of the ESPOIR cohort in 
which patients with RRP during the first year in the 
cohort, with a definition slightly different from the pre- 
vious one, showed increased structural damage progres- 
sion during the second and third years in the cohort [28]. 

The poor prognostic value associated with RRP is 
important and may be addressed by the development of 
prediction matrices to quantify the risk of RRP at one 
year in terms of baseline data. Such matrices are a tool 
than can identify patients with early RA who will show 
inadequate response to MTX or its equivalent [24,29-31]. 
These matrices have been developed in randomized clini- 
cal trials of patients with aggressive disease, who are not 
generally representative of patients with early RA. There- 
fore, the use of these matrices is limited in clinical prac- 
tice, as was recently shown in the BRASS cohort [32] . 

Thus, we conducted a study of data for the ESPOIR 
cohort, which enrolled patients with early arthritis from 
the community (with or without unfavorable prognostic 
factors). We aimed to develop a prediction algorithm and 
matrix to identify patients with early RA at risk of RRP 
despite early synthetic DMARD initiation. 

Materials and methods 

Patients 

Between December 2002 and March 2005, 813 patients 
with possible RA who were referred by rheumatologists 
and general practitioners to one of 14 regional centers 



were included in the ESPOIR cohort [33]. Inclusion cri- 
teria were age 18 to 70 years, more than two swollen joints 
for > 6 weeks and < 6 months, suspected or confirmed 
diagnosis of RA, and not taking any DMARDs or steroids 
except for < 2 weeks before enrollment. During the first 
year, patients were followed every six months. Clinical and 
biological data collected were disease activity by the Dis- 
ease Activity Score in 28 joints-4 variables (DAS28(ESR)- 
4v) [34] and functional ability by the Health Assessment 
Questionnaire (HAQ) [35]. Radiographs of the hand and 
feet (antero-posterior views) were obtained at each time, 
as was information on therapeutic regimen. Treatment 
strategies were not protocol-based in the ESPOIR cohort, 
and patients received usual care by their rheumatologist. 
The protocol of the ESPOIR Cohort study was approved 
by the ethics committee of Montpellier University Hospi- 
tal, France. All patients gave their signed informed consent 
to participate in the cohort. 

The current study involved data for ESPOIR patients 
with an RA diagnosis according to their rheumatologist 
and initiation of a first synthetic DMARD such as MTX 
or leflunomide with demonstrated structural efficacy for 
at least three months during the first year of follow-up in 
the cohort. 

Structural damage assessment 

X-ray data were collected in the radiography coordinating 
center and were read pair- wise by a well-trained investiga- 
tor blinded to clinical evaluation (intra-class correlation 
coefficient 0.99, SDD 0.966) [4]. Structural damage was 
assessed qualitatively by the presence of typical RA ero- 
sions, based on their location and aspect, and quantified 
according to the vSHS [36,37] . RRP was defined as change 
in vSHS (AvSHS) > 5 at 1 year [24,27,30,38]. 

Data management and statistical analysis 

Quantitative variables are expressed as mean ± SD and 
median. Qualitative data were expressed as number (%). 
The predictors of AvSHS were selected by a conventional 
two-step procedure. 

The Mann-Whitney U test (for numerical data) and 
Fisher's exact test (for categorical data) were used in a uni- 
variate analysis to establish the statistical significance of 
the relation between candidate predictors and RRP, with- 
out any a priori assumptions about the distribution of the 
variables. All variables with P < 0.1 were selected for mul- 
tivariate regression analysis. To construct the prediction 
matrix, quantitative variables selected in univariate analysis 
were categorized, the optimal threshold being selected on 
the basis of the variable distribution. 

The multivariate analysis was based on a logistic regres- 
sion model with a conventional backward stepwise proce- 
dure validated by a forward stepwise procedure whereby 
variables were optimized by the Akaike information 



Fautrel et al. Arthritis Research & Therapy 2012, 14:R249 
http://arthritis-research.eom/content/14/6/R249 



Page 3 of 9 



criteria, with P < 0.05. The overall discrimination power of 
the model was evaluated by receiver operating characteris- 
tic curve (ROC) analysis and the calculation of the area 
under the ROC (AUC). The fit of the model was assessed 
by the Hosmer-Lemeshow test. 

Matrix elaboration 

The RRP prediction matrix was developed by use of the 
model equation estimating the probability for one patient 
with early RA to display RRP at one year for each combi- 
nation of identified predictors. The 95% confidence inter- 
vals of individual probabilities were calculated with 1,000 
bootstrap replications, after removal of the outliers. 

All tests involved use of R 2.12.1 (R Foundation, 
Vienna, Austria). 

Results 

Patient characteristics 

From the ESPOIR cohort, 370 patients (45.5%) had started 
MTX [n = 335, mean dose 17.5 mg/week) or leflunomide 
(« = 35, mean dose 20 mg/week). These patients, referred 
to as 'synthetic DMARD (SD)-treated patients', were com- 
parable to the rest of the cohort in terms of age, sex, 



swollen and tender joints, and functional limitation. How- 
ever, they had shorter disease duration, higher biological 
inflammation values and more frequently were positive for 
rheumatoid factor (RF) or anti-citrullinated peptide anti- 
bodies (ACPA), had joint damage and satisfied the 2010 
criteria of the American College of Rheumatology/Eur- 
opean League Against Rheumatology (EULAR; Table 1). 
This indicates that the selected patients were more likely 
to have RA than those who did not receive one of the two 
drugs, illustrating a channeling bias. 

Patients who initially received MTX started the drug on 
average 27.4 ± 15.4 weeks after disease onset. Among 
them, 302 were still receiving MTX at 1 year, 19 had 
switched to another SD and 11 to a biologic agent, mainly 
TNF blockers (Table 2). Among patients receiving lefluno- 
mide, only two switched to a biologic agent at one year. 

Structural damage 

For the SD-treated patients, the mean structural progres- 
sion within the first year in terms of vSHS was 1.6 ± 5.5 
(median 0), mainly because of progression of joint erosion 
(Table 3, Additional file 1-Figure SI). The progression 
occurred in 126 (34.1%) patients (AvSHS > SDD) with 



Table 1 Baseline characteristics of the ESPOIR patients who received methotrexate or its equivalent as a first-line 
biologic agent for rheumatoid arthritis. 



SD-treated patients 
(number = 370) 



All ESPOIR patients 
(number = 813) 



Age, years 
Female sex 

Disease duration, weeks 
Swollen joint count in 28 joints 
Tender joint count in 28 joints 
ESR, mm/1 st hour 
CRP, mg/L 
DAS28(ESR)-4v 
IgM RF positivity 
ACPA positivity 

Typical erosion on radiographs 
vSHS score 

ACR/EULAR 2010 criteria 
HAQ score 
First-line agent: 

- No DMARDs 

- DMARDs without structural effect 

- Methotrexate or leflunomide 

- Other DMARDs with structural effect 

- Tumor necrosis factor blockers alone or in combination 



49.4 ±11.4 (51.5) 
271 (73.2%) 
15.2 ± 15.4 (57.7) 
7.9 ± 5.4 (7) 
8.7 ± 6.9 (7) 

32.7 ± 25 (26) 

24.8 ± 37.7 (11) 
5.4 ± 1.2 (5.2) 
204 (55.1%) 
185 (50%) 

66 (17.8%) 

6.02 ± 9.7 [54] 
316 (85.4%) 

1.03 ± 0.7 (1) 

n.a. 
n.a. 

370 (100%) 

n.a. 

n.a. 



48.1 ± 12.5 (50.1) 
624 (76.7%) 
31.6 ± 37.1 (21.3) 
7.2 ± 5.4 (6) 
8.4 ± 7 (6) 
29.4 ± 24.6 (22) 
20.3 ± 32.4 (9) 
5.1 ± 1.3 (5.1) 
376 (45.8%) 
315 (38.8%) 
1 00 (1 3.6%) 
3.71 ± 5.71 [54] 
582 (79.1%) 
0.979 ± 0.684 

207 (25.5%) 
117 (14.4%) 
396 (48.7%) 
56 (6.9%) 
37 (4.6%) 



Data are mean ± SD (median) or number (%); n.a.; not available. Baseline CRP level (normally < 10 mg/l), IgM and IgA RF (ELISA, Menarini, France; positive > 9 III/ 
ml) and anti-CCP2 antibodies (ACPA; ELISA, DiaSorin, France; positive > 50 U/ml) were detected in all patients with the same technique in a central lab (Paris- 
Bichat). ACPA, anti-citrullinated protein antibody; ACR/EULAR, American College of Rheumatology/European League Against Rheumatism; CRP, C-reactive protein; 
DAS28(ESR)-4v, Disease Activity Score in 28 joints-4 variables, using erythrocyte sedimentation rate; DMARD, disease-modifying anti-rheumatic drugs; ESR, 
erythrocyte sedimentation rate; HAQ, Health Assessment Questionnaire; IgM, immunoglobulin M; RF, rheumatoid factor; SD, synthetic DMARD; vSHS, van der 
Heijde-modified Sharp score. 
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Table 2 Evolution of DMARD treatment and disease during the first year for ESPOIR patients with early rheumatoid 
arthritis (RA). 



First DMARD 






Treatment at 1 year 








Type 


N 


Delay between RA onset and DMARD start 

(weeks) 3 




N 


DAS28(ESR)-4v at 
baseline 


DAS28(ESR)-4v at 1 
year 


Methotrexate 15 


335 


26.7 ±11.6 (24.3) 


Methotrexate 


302 


5.3 ± 1.3 (5.2) 


3.2 ± 2.7 (3.0) 






35 ± 15.3 (39.4) 


Other synthetic 
DMARD(s) c 


19 


5.3 ± 1.4 (5.3) 


4.1 ± 1.5 (4.1) 






30.5 ± 0.7 (21.6) 


Biologic agent d 


11 


5.7 ± 0.4 (5.5) 


3.9 ± 0.9 (4.2) 


Leflunomide 


35 


48.9 ± 1 1 .6 (46.9) 


Methotrexate 


9 


5.8 ± 1 .3 (5.6) 


4.0 ± 2.7 (2.9) 






22.5 ± 15.3 (20.7) 


Other synthetic 
DMARD(s) 


24 


5.4 ± 1 .4 (5.5) 


3.3 ± 1.5 (3.0) 






37.8 ± 0.7 (37.8) 


Biologic agent d 


2 


5.9 ± 0.4 (5.9) 


4.2 ± 0.9 (4.2) 



Data are mean ± SD (median). a Delay between RA onset and DMARD start (weeks) (P = 0.479): Methotrexate: 27.2 ± 15.1 (24.7); Leflunomide: 30.1 ± 18.1 (24.6); 
b at one year: one patient did not receive any DMARDs and data for two patients were not available; c leflunomide or salazopyrine; d adalinumab, etanercept, 
infliximab or anakinra. DAS28(ESR)-4v, Disease Activity Score in 28 joints-4 variables, using erythrocyte sedimentation rate; DMARD, disease-modifying anti- 
rheumatic drugs; N, number. 



mean vSHS change of 5.4 ± 7.0 (median 2). Among them, 
41 had RRP, which represented 11.1% of all SD-treated 
patients and 32.2% of patients with disease progression. 

RRP determinants 

The variables associated with RRP on univariate analysis 
were disease duration, swollen joint count (SJC) and CRP 
level (only when considered a categorical variable), RF and 
ACPA status, alone or combined, and at least one typical 
erosion seen on hand or foot radiographs (Table 4). Of 
note, there were no differences between the two treatment 
populations (MTX or leflunomide) with regard to these 
determinants. The final multivariate model included only 
SJC and CRP level (with two different thresholds for each), 
as well as combined RF or ACPA status and presence of 
typical RA erosions, which were predictors of RRP 
(Table 5). ROC analysis confirmed the good discriminating 
power of the model with an area under the curve (AOC) 
of 0.754 [see Additional file 1, Figure S2]. 

Matrix elaboration 

From these data, we developed a risk matrix (Figure 1, 
Additional file 1, Figure S3) that indicates the probability 
of RRP for patient profiles in terms of baseline SJC, CRP 



level, ACPA status, and erosions seen on radiography. 
For example, for a patient without typical RA erosions on 
radiographs, no ACPA positivity, who receives MTX, the 
probability of RRP at one year is 2% with SJC < 14 and 
CRP level < 4 mg/L but 14% with SJC > 20 and CRP level 

> 35 mg/L. However, for a patient with typical RA ero- 
sions on radiographs, ACPA positivity, who receives 
MTX, the probability of RRP is 18% with SJC < 14 and 
CRP level < 4 mg/L but 64% with SJC > 20 and CRP level 

> 35 mg/L. In terms of the high RRP probability, treat- 
ment for the latter patients with only MTX or equivalent 
could be considered suboptimal. Such patients repre- 
sented 0.5% (n = 2) of the ESPOIR population if the 50%- 
threshold is used (red boxes) and up to 12.4% (n = 46) if 
the 25%-threshold is used (red and orange boxes). 

Discussion 

The present study allowed for the construction of a 
matrix to predict the risk of RRP for patients with early 
RA despite MTX or leflunomide therapy. Since RA is 
considered a medical emergency [39] requiring rapid 
referral to a rheumatologist [1-3] and early initiation of 
DMARDs to prevent disease progression, the ESPOIR 
matrix might help rheumatologists in daily practice 



Table 3 Structural disease progression in patients with early RA who received DMARDs and the entire ESPOIR cohort. 



SD-treated patients All ESPOIR patients 



All patients 


(number = 370) 


(number = 736 a ) 


- AvSHS total 


1 .6 ± 5.5 (0) 


0.96 ± 4.4 (0) 


- AvSHS erosion 


1 .5 ± 4.7 (0) 


0.96 ± 3.8 (0) 


- AvSHS narrowing 


0.2 ± 1 .8 (0) 


0.01 ± 1.81 (0) 


Patients with structural progression b 


(number = 126, 34.1%) 


(number = 187, 25.4%) 


-AvSHS total 


5.4 ± 7.0 [54] 


5.0 ± 6.3 [54] 


- AvSHS > 5 points (RRP) 


41 (11.1%) 


58 (7.9%) 



Data are mean ± SD (median) or number (%). a X-ray data were available for only 736 patients; Structural progression was defined as AvSHS > 1 point (that is, the 
smallest detectable difference). DMARD, disease-modifying anti-rheumatic drugs; RRP, rapid radiographic progression; SD, synthetic DMARD; AvSHS, change in 
van der Heijde-modified Sharp score. 
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Table 4 Association of main baseline characteristics of patients with early RA and rapid radiographic progression 



(RRP) of rheumatoid arthritis (univariate analysis). 





With RRP 
(number = 41) 


Without RRP 
(number = 329) 


P value 3 


Age, years 


49.8 ± 1 2 (52.7) 


49.3 ± 1 1 .5 (7) 


0.6 


Sex 


31 (75.6%) 


240 (73%) 


0.85 


Disease duration, weeks 


18.6 ± 8.4 (18.6) 


14.7 ± 8.1 (13) 


0.007 


SJC in 28 joints 


8.6 ± 6.1 (8) 


7.8 ± 5.3 (7) 


0.5 


< 14 


31 (75.6) 


285 (86.6) 




14 to 20 


6 (14.6) 


32 (9.7) 


0.08 


> 20 


4 (12.2) 


1 2 (3.6) 




TJC in 28 joints 


8.3 ± 6.3 (7) 


8.7 ± 6.9 (7) 


0.97 


ESR, mm/1 hr 


32.6 ±21.3 (30) 


32.7 ± 25.4 (25) 


0.5 


CRP, mg/L 


26.2 ± 27.9 (14) 


24.6 ± 38.7 (11) 


0.12 


< 4 


4 (9.7) 


89 (27.1) 




4 to 35 


25 (61) 


168 (51.1) 


0.04 


> 35 


12 (29.3.8) 


72 (21.9) 




Elevated ESR or CRP level 


38 (90.3) 


289 (87.8) 


0.45 


DAS28(ESR)-4v 


5.3 ± 1 .2 (5.4) 


5.3 ± 1.2 (5.2) 


0.8 


RF positivity 


29 (70.7) 


1 75 (53.2) 


0.04 


ACPA positivity 


31 (75.6) 


1 54 (46.8) 


0.0008 


RF or ACPA positivity 


32 (78) 


190 (57.8) 


0.01 


HAQ score 


0.95 ± 0.6 (1) 


1.04 ± 0.7 (1) 


0.5 


Typical RA erosion 


1 8 (44) 


48 (14.6) 


< 0.0001 


Prednisone > 7.5 mg/d 


2 (4.9) 


35 (10.6) 


0.4 


> 5 mg/d 


7 (17.7) 


91 (27.7) 


0.19 


Delay before 1st DMARD initiation > 6 months after RA onset 


22 (53.7) 


147 (44.7) 


0.32 


Satisfaction of 2010 ACR/EULAR criteria 


38 (92.7) 


278 (84.5) 


0.24 



Data are mean ± SD (median) or number (%). RRP-positive status was defined by progression of the vSHS total score > 5 points between baseline and one year 
(of note, among the 46 patients with RRP, 44 had progression with vSHS erosion score > 5 points and 13 with vSHS narrowing score > 5 points). a Student f test 
or Mann-Whitney U test for continuous variables and Fischer exact test for categorical variables, with P < 0.1 . ACPA, anti-citrullinated protein antibody; ACR/ 
EULAR, American College of Rheumatology/European League Against Rheumatism; CRP, C-reactive protein; DAS28(ESR)-4v, Disease Activity Score in 28 joints-4 
variables, using erythrocyte sedimentation rate; DMARD, disease-modifying anti-rheumatic drugs; ESR, erythrocyte sedimentation rate; HAQ, Health Assessment 
Questionnaire; RF, rheumatoid factor; SJC, swollen joint count; TJC, tender joint count; vSHS, van der Heijde-modified Sharp score. 



identify patients at high risk of SD failure, and thus 
make better-informed and evidence-based therapeutic 
decisions. The community-based nature of the study is 
one of the major strengths of this work, which makes 



Table 5 Determinants of rapid radiographic progression 
of RA for patients with early RA who received DMARDs 
(multivariate analysis). 





Estimate 


Standard 
Error 


z 

value 


P 

value 


Swollen joint count 14 to 
20 


0.27 


0.53 


0.52 


0.60 


Swollen joint count > 20 


1.25 


0.68 


1.84 


0.06 


CRP 4 to 35 mg/L 


0.83 


0.57 


1.45 


0.15 


CRP > 35 mg/L 


0.86 


0.63 


1.36 


0.17 


ACPA status 


1.11 


0.40 


2.75 


0.006 


Typical RA erosion 


1.31 


0.63 


3.53 


0.0004 


Intercept 


-3.94 


0.58 


-6.74 


1 .62e- 
11 



ACPA, anti-citrullinated protein antibody; CRP, C-reactive protein; DMARD, 
disease-modifying anti-rheumatic drugs; RA, rheumatoid arthritis 



patient typology more representative of daily clinical 
practice than that of randomized clinical trials. The 
ESPOIR cohort included patients with early arthritis and 
possible, probable or definite diagnosis of RA. However, 
no specific disease activity level or prognostic markers 
of severity were required for inclusion [33] . 

For several reasons, we focused on patients who had 
started receiving MTX or leflunomide during the observa- 
tion period. First, therapeutic decisions were not protocol- 
based in the ESPOIR cohort, and patients received treat- 
ment according to standard care by their rheumatologist. 
Patients with the most active or structurally aggressive dis- 
ease were more likely to receive the most effective drugs 
such as MTX or TNF blockers, which highlights a chan- 
neling bias [4,40] . Four gross therapeutic strategies can be 
distinguished: no DMARD during the first year in the 
cohort, DMARDs without demonstrated structural efficacy 
such as hydroxychloroquine, DMARDs with demonstrated 
structural benefit such as MTX or leflunomide, or biolo- 
gies such as TNF blockers. Without any adjustment, the 
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Absence of typical RA erosion on radiographs 



Presence of typical RA erosions on radiographs 



AC PA 
positivity 



AC PA 
negativity 



CRP > 35 



4 < CRP < 35 



CRP<4 



CRP > 35 



4 < CRP < 35 



CRP<4 



SJC<14 


14 < SJC < 20 


SJC > 20 


SJC < 14 


14 < SJC < 20 


SJC > 20 


0.12 


0.16 


0.33 


0.34 


0.40 




[0.05; 0.23] 


[0;0.91] 


[0.10; 0.84] 


[0.17; 0.54] 


[0.16; 0.67] 


[0 28 6 0 4 88] 


0.12 


0.16 


0.32 


0.33 


0.40 


0.64 


[0.07; 0.19] 


[0;0.9] 


[0.09; 0.63] 


[0.19; 0.50] 


[0.18; 0.52] 


[0.24; 0.88] 


0.06 


0.07 


0.17 


0.18 


0.22 


0.43 


[0.01; 0.13] 


[0.01 ; 0.24] 


[0.03; 0.52] 


[0.04; 0.35] 


[0.04; 0.52] 


[0.07; 0.82] 


0.04 


0.06 


0.14 


0.15 


0.18 


0.37 


[0.01; 0.10] 


[0.01; 0.16] 


[0.03; 0.33] 


[0.04; 0.31] 


[0.05; 0.43] 


[0.08; 0.68] 


0.04 


0.06 


0.14 


0.14 


0.18 


0.37 


[0.01 ; 0.08] 


[0.01; 0.15] 


[0.03; 0.33] 


[0.04; 0.28] 


[0.05; 0.44] 


[0.08; 0.69] 


0.02 


0.03 


0.07 


0.07 


0.08 


0.20 


[0; 0.04] 


[0; 0.08] 


[0.01; 0.19] 


[0.02; 0.15] 


[0.01 ; 0.24] 


[0.02; 0.52] 



I 



RRP Risk 
a 50% 
25 < <50% 
10 < <25% 
< 10% 



Figure 1 ESPOIR prediction matrix. ESPOIR prediction matrix for use in daily practice in assessing the risk of rapid radiographic progression 
(RRP; change in vSHS > 5 points at 1 year) in patients with early rheumatoid arthritis in terms of baseline characteristics. ACPA, anti-citrullinated 
protein antibody; CRP, C-reactive protein (mg/L); vSHS, van der Heijde-modified Sharp score. 



progression in vSHS was substantially higher in the latter 
two groups, which shows that patients receiving the most 
efficacious drugs were identified by their rheumatologists 
as being likely to have active and erosive RA [4] . The use 
of propensity scores could in theory adjust on these con- 
founding factors. Although different models were tested, 
we did not achieve a satisfactory adjustment, able to make 
the different treatment groups comparable at baseline. 
Thus, the analysis was stratified and only focused on the 
more clinically relevant and adequately sample-sized 
group, that is, the group of patients treated with either 
MTX or leflunomide, the two most widely used therapies 
in early RA. Although several synthetic DMARDs have 
shown their ability to prevent structural damage progres- 
sion, only two seem to have similar efficacy (that is, MTX 
and leflunomide) in a recent systematic literature review 
conducted under the auspices of EULAR [15]. Therefore, 
sulfasalazine, gold salts or other synthetic DMARDs were 
excluded in the patient selection process. This choice also 
appeared more clinically relevant because this therapeutic 
option is considered optimal for most patients [17]. Our 
data indicate that other options, for example, biologies, 
could be better for the patients with high RRP risk (red 
boxes in Figure 1). In an explanatory analysis, the patients 
treated with TNF-blockers as first line agent {n = 37, with 
complete data available for only 27) were plotted in the 
matrix; their observed RRP risk was lower than the 
expected RRP risk. 

Several prediction matrices have been recently proposed. 
Three were derived from results of randomized controlled 
clinical trials assessing the efficacy of TNF-blocking agents 
[20,24,29-31,41], so patients with early RA with substantial 
disease activity and/ or severity were included as candidates 
for biologies. Moreover, in these trials, therapeutic adapta- 
tion was protocol-based - fixed strategy over one year or 
tight control according to a predefined scheme - which, in 



both cases, is not quite consistent with usual care in which 
treatment adaptations are often looser, thus enabling sub- 
optimal disease control [42,43] . Only one algorithm was 
developed from data from an early RA cohort, SONORA, 
in Canada and the United States but has not been pub- 
lished [44] . The predictors identified to build the different 
matrices were partly overlapping. In the ASPIRE matrix 
[24], the three constitutive variables comprising the matrix 
were RF status, SJC and erythrocyte sedimentation rate 
(ESR); baseline joint damage was not included, although it 
is a strong predictor of further structural damage in RA 
[45]. In the BeSt study [30], RF or ACPA status, number 
of erosions at baseline and CRP level but not SJC, another 
marker of disease severity, were used in the risk matrix 
[45]. In the SWEFOT trial [20,29,31,41], the proposed 
matrix included smoking status, CRP level and erosions at 
baseline; ACPA status was not significantly associated 
with RRP, which is not consistent with the literature [45]. 
This finding may be explained in part by the confounding 
of the strong association of smoking status and ACPA sta- 
tus [46-48]. The ESPOIR matrix may appear more com- 
prehensive than the other matrices because it included all 
the known RA prognostic markers. Of interest, in the 
observational cohort SONORA, the proposed matrix 
included most of the same parameters: baseline DAS28 
(ESR)-4v , vSHS and ACPA status [44]. 

Our study has some limitations. The first pertains to 
MTX dosage, because several doses have been proposed in 
the literature. Our patients received 17.5 mg/week, on 
average, which may be considered slightly lower than that 
used in other settings [49]. Another concern relates to 
low-dose prednisone therapy, which has been efficacious 
in preventing structural damage progression [50-52]. In 
the ESPOIR cohort, 42% of the patients received steroids 
during the first year of follow-up, with a mean dose of 5.5 
± 3.7 mg/d (median 5 mg/d). This daily dose was quite 
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low, and our analyses did not reveal any association 
between prednisone intake and RRP. 

The quality of care could not be assessed in the 
ESPOIR cohort, especially in terms of optimal or subopti- 
mal timing of RA diagnosis or MTX introduction, two 
key elements of future disease control [3,4,13,53]. More- 
over, the implementation and respect of tight control was 
not identifiable in the cohort and may have had an 
impact on structural progression [42]. However, this 
situation may suggest that our results are representative 
of clinical practice, as stated above. Another important 
limitation is the lack of validation of our matrix in a dif- 
ferent population, which is a limitation in all RA 
matrices. Attempts have been made to cross-validate the 
ASPIRE, BeSt and SWEFOT matrices in the SWEFOT 
trial population [41] and in the BRASS cohort, an estab- 
lished RA cohort [32]. In both cases, the performance of 
the different matrices in different populations was disap- 
pointing. Whether the more comprehensive ESPOIR 
matrix may be optimal in different populations remains 
to be demonstrated, as well as the performance of the 
other matrices in the ESPOIR population (work in pro- 
gress). Finally, all matrices have limitations inherent in 
their development that relate to multivariate analyses 
performed in limited samples of patients. The RRP risk 
was thus a predicted risk as opposed to an observed risk. 
Some of the matrix boxes were rather poorly populated 
and may explain the only fair performance of matrices in 
populations different from those in which they have been 
developed. 

Conclusions 

The ESPOIR matrix is a novel tool developed in a real life 
setting to help rheumatologists in usual daily practice to 
identify patients with early RA at high risk of RRP despite 
MTX or leflunomide therapy. The performance of this 
tool needs to be validated in other patient populations. 
Then, it will become a quite relevant instrument to guide 
rheumatologists in their therapeutic decision making, 
especially to detect patients for whom MTX may not be 
optimal therapy. 

Additional material 
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